Sign-Alternating Interaction Mediated by Strongly-Correlated Lattice Bosons 
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We reveal a generic mechanism of generating sign-alternating inter-site interactions mediated by 
strongly correlated lattice bosons. The ground state phase diagram of the two-component hard-core 
Bose-Hubbard model on a square lattice at half- integer filling factor for each component, obtained 
by worm algorithm Monte Carlo simulations, is strongly modified by these interactions and features 
the solid-|-superfluid phase for strong anisotropy between the hopping amplitudes. The new phase 
is a direct consequence of the effective nearest-neighbor repulsion between "heavy" atoms mediated 
by the "light" superfluid component. Due to their sign-alternating character, mediated interactions 
lead to a rich variety of yet to be discovered quantum phases. 
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The first proposal for studying models of strongly cor- 
related systems with cold atoms in optical lattices was 
put forward a decade ago Since then, control over 
lattice geometry and interaction strength has increased 
dramatically, opening up new directions in the studyof 
quantum phases of cold gases. (For reviews, see 
Thanks to refinements in experimental and theoretical 
tools, it is now possible to look at exotic quantum states 
which arise in bosonic systems with pseudospin degrees 
of freedom or multiple species. In the realm of two- 
component systems, one goal is to realize models of quan- 
tum magnetism by using hyperfine states of an atom 
By controlling superexchange interactions of particles 
confined in an optical lattice, it is possible to switch 
between different ground states Another important 
development is experimental realization of heteronuclear 
bosonic mixtures of ^''Rb — '^^K in a three-dimensional 
optical lattice 6]. Moreover, Ref. reports results for 
fine control over interspecies scattering length, includ- 
ing the zero-crossing point. These achievements indi- 
cate that two-component systems in optical lattices with 
tunable interspecies interaction via Feshbach resonances 
are within the reach of current experiments. The two- 
component 2D bosonic system is also in the focus of Op- 
tical Lattice Emulator project supported by DARPA and 
aimed at the development, within the next few years, 
of experimental tools of accurately mapping phase dia- 
grams of lattice systems by emulating them with ultra- 
cold atoms in optical lattices. 

Experimental studies of lattice solids — states with bro- 
ken translation symmetry — are intriguing and fundamen- 
tally important, especially in 2D. For the prominent ex- 
ample, we refer to the problem of deconfined critical- 
ity proposed for the solid-to-superfluid quantum phase 
transition in 2D. The matter of interdisciplinary interest 
here is to validate the idea of a (hidden) duality between 
the superfluid and solid orders leading to a conceptually 
new criticality ^SJ. Lattice solids also offer the possibil- 
ity of having a supersolid phase featuring both broken 
translation symmetry and the ability to support a super- 



flow, e.g. in a single-species square-lattice bosonic sys- 
tem with soft-core on-site interactions^ and appropriately 
strong nearest-neighbor interactions [l3j . 

Solid phases in the single-species bosonic system re- 
quire going beyond the on-site interaction. The standard 
Bose-Hubbard model ^ supports only two phases: a su- 
perfluid (SF) and a Mott insulator (MI); the latter is 
not a solid because it lacks the broken translation sym- 
metry. A considerable theoretical effort has been made 
to understand how inter-site interactions can be gener- 
ated in atomic gases. One proposal is to use cold polar 
molecules [l^ featuring long range dipole-dipole interac- 
tions (se e also review and references therein). In 
Refs. [la], the authors suggest a technique for tuning the 
shape of long-range interactions between polar molecules 
by applying static and microwave fields. Another, exper- 
imentally more challenging, proposal is to excite atoms 
to higher bands in an optical lattice [iTj. 

Within this framework, the two-component bosonic 
system with purely on-site interactions is a reasonable 
alternative route to obtain exotic single-species systems 
(at present, it is hard to reach low temperatures with 
lattice fermions). At a commensurate filling and strong 
enough interaction, opening a gap in the net-charge sec- 
tor, the two-component mixture becomes equivalent to a 
single-component system with nearest-neighbor interac- 
tions, describing the iso-spin sector 9, 10, 11, 12). The 
checker-board (CB) solid phase arising in this case is 
equivalent to the Neel antiferromagnet [§, [13, [i2| • The 
CB-SF quantum phase transition (with respect to the 
original components, the SF state is a super-counter-fluid 
(SCF) [lo| . or, equivalently, planar ferromagnet in iso- 
spin terminology [l^, [l3|) is known to be of the first or- 
der. Unfortunately, the supersolid phase is not predicted 
in this parameter regime. There have been extensive 
theoretical studies of quantum phases in two-component 
systems [13, [Hi, [H, 111 [III, but, to the best of our 
knowledge, all studies overlooked the possibility of having 
various solids in the absence of the net-charge localiza- 
tion. Namely, they missed a generic mechanism of indue- 
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ing inter-site sign-alternating interactions by strongly- 
correlated bosonic environment (cf. (20|'). analogous to 
the RKKY (Ruderman-Kittel-Kasuya-Yosida) interac- 
tion mediated by fermions f2l|. 

In this Letter, we reveal and quantify the mechanism 
of the mediated sign-alternating interactions, and discuss 
it in the context of the ground state phase diagram for 
hard-core bosons with repulsive interspecies interaction 
at half filling for each component: 

^ = - E (^aal a, + t^b] 6, ) + C/ E ) , (1) 

<ij> i 

where, a|(ai) and are bosonic creation (annihi- 

lation) operators, ia and th are hopping matrix ele- 
ments between the nearest neighbor sites for two species 
of bosons {A and B) on a simple square lattice with 

N = LxL sites, and nl'^^ = ajai, nf"^ = bjbi. This model 
can be implemented experimentally and is considered 
to be the simplest one with purely contact interactions 
and yet highly nontrivial phase diagram. 

Model ^ was studied previously using a combination 
of variational and mean field theories [I4I which, in gen- 
eral, can not guarantee the accuracy of results. With 
Monte Carlo (MC) simulations by Worm Algorithm , 
we obtain the first precise data for the ground state phase 
diagram. For weak anisotropy between ia and t^, and 
large U, our results confirm the basic phases and transi- 
tions between them proposed in Ref. [ij]- We, however, 
find strong quantitative differences (up to 50% to 100%) 
in the location of transition lines. For large anisotropy 
and moderate-to-weak interactions we find a completely 
new structure of the phase diagram. It is shaped by 
the effective Hamiltonian obtained for the "heavy" (small 
hopping) component after the "light" component is inte- 
grated out. The resulting nearest-neighbor and longer- 
range interactions (similar to the effective potential be- 
tween the ions in solids mediated by electrons) stabilize 
the checker-board (CB) solid phase of heavy atoms for 
sufficiently strong anisotropy between and t^. A sur- 
prising result of the present study is that effective me- 
diated interactions are oscillating from strong on-site at- 
traction to much weaker nearest neighbor repulsion and 
back to a tiny attractive tail. In a broad perspective, 
this type of mediated interactions will result in interest- 
ing solid, and guaranteed supersolid [T3| orders in related 
models. Moreover, for soft-core bosons, one can look for 
phases and phase transitions which involve multi-particle 
bound states and order parameters ( "multi-mers" ) . 

Before we discuss our findings in more detail let us 
review the key phases and limiting cases of model ([1]). 
In the strong coupling limit, U ^ ta_,t\^, it can be 
mapped (within the second-order perturbation theory) 
onto the spin- 1/2 Hamiltonian (see e.g. El [11) 
Hxxz = T,<ij>[-Jxy{c^j<^f+'^j<^!) + Jz'^j'^!] with pos- 
itive Jxy, Jz ^ /U. The latter features two possible 
ground states: (i) an antiferromagnetic state with z- 
Neel order for Jz > Jxy, and (i) an Xl"- ferromagnetic 



state for J^y > Jz- In bosonic language, the z-Neel 
state corresponds to the CB solid order for both A- 
and B-particles (we will abbreviate it as 2CB). It is 
characterized by non-zero structure factor S'a,b(k) — 

iV~^ exp [ikr] (no'^'^^rir*^'^^). The XF-state is repre- 
senting the SCF phase featuring an order parameter 
(at6). Both 2CB and SCF have to be regarded as Mott 
insulators as far as the total number of particles is con- 
cerned, i.e. there exist a finite gap to dope the system. 
Thus only counter-propagating A- and B-currents with 
the zero net particle flux posses supcrfiuid properties 
in the SCF state. Under the mapping one finds that 
Jz > Jxy everywhere except at ta = when the spin 
Hamiltonian becomes SU(2)-symmetric. Thus higher- 
order symmetry-breaking terms are necessaxy to decide 
which phase, 2CB or SCF, survives. Ref. [iJl provided 
a variational argument showing that SCF is stabilized in 
the vicinity of the <a = line, and our data unambigu- 
ously confirm the validity of this conclusion. 

At weak inter-species interaction, U ^ ta,t\,, we 
expect that the ground state is that of two miscible 
strongly interacting (due to hard-core intra-component 
repulsion) superfluids (2SF). Finally, for % <^ U < t^, 
we should have a phase where B-particles form the CB 
solid if effective interactions mediated by the superfluid 
A-component are repulsive and short-ranged (we abbre- 
viate it as CB-HSF). 

Our simulation method is based on the lattice path 
integral representation and Worm Algorithm [22|. The 
original version was generalized to deal with two- 
component systems following ideas introduced for clas- 
sical j-current models (23| . The simulation configura- 
tion space now includes the possibility of having two 
types of disconnected worldlines (worms) representing 
off-diagonal correlation functions (Green's functions). In 
order to allow efficient sampling of the SCF phase (any 
paired phases for that matter) it is necessary to enlarge 
the configuration space and consider worldline trajecto- 
ries with two worms propagating simultaneously. The 
results for the phase diagram are summarized in Fig. [T] 
To detect the SCF phase we have calculated the stiffness 
of the relative superfluid flow from the standard winding 
number formula psCF = /3"^ ((Wa — Wb)^), where 
Wa(b) are winding numbers of worldlines A{B), and /3 
is the inverse temperature. In SCF the sum of winding 
numbers is zero in the thermodynamic limit. We confirm 
the SCF ground state for ta ^ and sufficiently strong 
interactions. It survives at arbitrary large U along the 
diagonal ia = ^b, directly demonstrating that higher or- 
der terms in the effective spin- 1/2 Hamiltonian break the 
SU(2) symmetry in favor of the XY-order. To locate the 
weakly first-order superfluid-solid 2CB-SCF line (circles 
in Fig.[T]) we have used the flowgram method which works 
well for both first and second order transitions, and is 
particularly helpful for telling the former from the latter 
(see Ref. [23] for details). 

Though the 2CB-SCF transition is expected to be first- 
order, one may not exclude the possibility of the interme- 
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FIG. 1: (Color online). Phase diagram of model ([T| on a 
square lattice at half-integer filling factor for each compo- 
nent. The observed transition lines are 2CB-SCF (first-order), 
SCF-2SF (second-order), 2CB-2SF (first-order), 2CB-CB-HSF 
(second-order), and CB+SF-2SF (first-order). Lines are used 
to guide an eye. 




FIG. 2: The C(r) function in a system of light hard-core 
bosons at half-integer filling; the distance r is measured in 
units of lattice spacing. The calculation was done for the L x 
L = 10 X 10 system at low temperature. The nearest neighbor 
repulsion is clearly visible, though the overall strength of the 
effective coupling is very small. Error bars are smaller than 
symbol size. 



diate supersolid phase. We did search for evidence of the 
state featuring both the SCF and CB orders but did not 
find any. On the contrary, we observed phase coexistence 
which brings us to the conclusion that the 2CB-SCF line 
remains un-split. 

The continuous SCF-2SF transition (squares in Fig.[T]), 



is expected to be in the {d + l)-dimensional U(l) univer- 
sality class characteristic of the MI-SF quantum phase 
transition |23!i| . As the system crosses into the 2SF 
phase, it develops single-component order parameters 
(a) 7^ and (b) ^ along with the non-zero su- 
perfluid stiffness in the total winding number channel, 
P2SF = /3~^((Wa -t- Wb)^). To locate the transition line 
precisely we have employed standard finite size scaling 
arguments and extracted the critical point from the in- 
tersection of P2SF{U/ta^b)L curves calculated for different 
system sizes L at /3 oc L. 

The SCF phase disappears for U/ta^t, < 8. In this 
parameter region, the system undergoes the first-order 
2CB-2SF transition (stars in Fig. [T]) up to C//ia,b ^ 4 
where the 2CB phase disappears. 

Not too surprisingly for the 2D case, the mean-field 
and variational treatments turn out to rather inaccurate 
quantitatively: The actual transition lines are 50% to 
100% away from the predicted ones. For com par ison, in 
Fig. [T] we have used the same units as in Ref. (2 = 4 
is the coordination number on the square lattice). Even 
for strong inter-species interaction J7/ta,b ~ 16 we do not 
find good quantitative agreement. 

We now turn to the most interesting, for purposes 
of this Letter, region of the phase diagram with strong 
anisotropy between hopping amplitudes. In this regime 
one of the components is much heavier (let it be com- 
ponent B) than the other. As the hopping amplitude ia 
increases, the light component undergoes a second-order 
MI-SF transition in {d+ l)-dimensional U(l) universality 
class (diamonds in Fig. [1]). [Since translation invariance 
is broken by the CB order of B-particles, the filling factor 
of A-particles is unity per unit cell.] Beyond this transi- 
tion, the superfiuid A-bosons provide an effective interac- 
tion for B-bosons which for sufficiently small tb stabilizes 
the CB order in the heavy component. Though formally 
CB-I-SF breaks both translation and gauge symmetries it 
can not be called a supersolid because the density wave 
in the A-component is induced from "outside" by an in- 
sulating heavy phase (this is reminiscent of conventional 
solid-state superconductors). However, in cold atomic 
systems with A and B particles referring to different hy- 
perfine states of the same atom we have an experimental 
"knob" to render the CB-f SF phase a genuine supersolid 
by inducing an arbitrary small hybridizing interaction 
g{a^b + H.c.). 

To describe the CB-I-SF phase semi-quantitatively let 
us look at the limit ib ■C C/ ^ ta. Within the lowest- 
order perturbation theory in U an effective mediated in- 
teraction between heavy atoms is pairwise and of the or- 
der of U'^ / zta,. Since the superfiuid liquid of hard-core 
A-particles is strongly correlated, the actual strength 
and form of Vfj^ has to be computed numerically. We 
determine it from the exact (at U/ta. — > 0) relation 
V^ir) = (C/Vfa)C(r), where 



C(r) 



lim 



n{r) - 1/2 
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is the scaled density profile in response to the heavy 
atom at the origin, and r is the dimensionless distance 
between the B-atoms. Since C{r) is expected to decay 
exponentially fast with r (26j . the simulated system can 
be relatively small. The result is shown in Fig. [2l It 
turns out that V°^{r) is more than an order of magni- 
tude smaller then a naive estimate / zta- Moreover, it 
is sign-alternating, with strong on-site attraction, much 
weaker nearest-neighbor repulsion, and nearly negligible 
attractive tail for longer range interactions. For the hard- 
core model studied here, the dominant interaction is the 
nearest-neighbor one, ^^{1) 5.4 x 10~^ C^^Aa, ex- 
plaining the origin of the CB phase of heavy atoms. 

On the basis of the melting first-order 

(CB-|-SF)-2SF transition is predicted to occur along the 
th « V''^{l)/2, or AzHht^/U^ « 0.17, line, i.e. the 
CB+SF phase survives close to the vertical and horizon- 
tal axis in Fig. [T] all the way to ta,b/C/ — > oo. The above 
asymptotic estimate is less than a factor of two smaller 
than the data points in Fig. [1] even though U is still 
relatively large. 

In conclusion, we have presented accurate results, 
based on path integral Monte Carlo, for the phase di- 
agram of the two component hard core Bose-Hubbard 
model on a square lattice and half-integer filling factor 
for each component. The system can be realized experi- 



mentally with heteronuclear bosonic mixtures in optical 
lattices with tunable interspecies interactions. We reveal 
the existence of an additional CB-I-SF state which radi- 
cally changed the topology of the phase diagram. The 
CB+SF phase, which exists for strong anisotropy be- 
tween the hopping amplitudes and weak enough inter- 
action, is a direct consequence of effective interactions 
mediated by the light, strongly correlated superfluid com- 
ponent. Mediated interactions are sign- alternating, and 
lead to exciting possibilities of realizing new quantum 
phases in the two-species bosonic systems. 

Many questions remain open. Finite temperature 
properties and melting of the z-Neel/a;?/-ferromagnet 
phases are of interest in the study of quantum magnetism. 
Studies of soft core bosons are of special interest because 
they admit the possibility of forming "multi-mers" due to 
strong mediated on-site attraction. One may also study 
supersolid phases on square and triangular lattices, not 
to mention the need for dealing with more realistic sys- 
tems, i.e. including effects of parabolic confinement and 
finite number of particles as in experimental setups. 
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